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1. Introduction 

The production of gravitational waves by binary systems, in the limit in which mass and 
momentum transfer do not occur, has been well described with high accurate precision 
with both the Post-Newtonian and the Post-Minkowskian approaches. These are the 
most useful tools in general relativity to provide good estimates of the features of the 
gravitational wave sources. Recently, great advances in these methods have improved 
their precision, using different multipole expansions. For a review of these issues we 
refer the reader to [1]. 

On the other hand, the characteristic formulation of general relativity on outgoing 
null cones provides an interesting point of view for treating gravitationally radiating 
systems, since it is based on the radiation’s coordinates. Originally formulated by Bondi 
et. al. in the 1960s pa e], it has served as the basis for the development of recent and 
accurate characteristic numerical codes. These have evolved from the quasi-spherical 
regime, where the non-linear angular terms are neglected [4j, to the full non-linear case 
[S]. With the inclusion of (perfect fluids) matter sources into the numerical codes [6], 
post hoc extensions were made, putting the system of equations in terms of only first 
order angular derivatives [7]. The development of the Cauchy-characteristic matching 
(CCM) method [8] allows the implementation of highly accurate gravitational wave 
propagation techniques. Great advances in the precision for the angular derivatives 
were made through the introduction of a gnomonic atlas to cover the angular manifold 
PED]. Also, there has been impressive progress in spectral methods m and in mixed 
methods, combining finite-difference schemes and spectral representation of the angular 
part of the field equations, including use of the Runge-Kutta 4 method to perform the 
time integration paiiEuii]. 

The main purpose of such codes, after all, is to simulate the gravitational radiation 
patterns produced, for example, in stellar collapses or by coalescing black hole—black 
hole (BH-BH), black hole—neutron star (BH-NS) or neutron star—neutron star (NS- 
NS) binary systems. As calibration tools, some toy models have been constructed 
from analytical solutions of the Einstein field equations in the linear regime (see, e.g., 

mmmm)- 

These toy models allow one to study the boundary conditions that the metric and their 
derivatives should satisfy in different situations in a simple, clear and straightforward 
way. Bishop et. al. Pi EH and Kubeka [18], provide a prescription to solve analytically 
the field equations for the vacuum in the linear regime. Then, they construct solutions 
for a variety of simple systems, such as a static (dynamical) thin shell, and a binary 
system composed of point particles of equal masses moving in circular orbits. However, 
the issue concerning the boundary conditions deserves additional studies, and most 
importantly, binaries of different masses in circular orbits have not been considered and 
therefore deserve closer scrutiny. 

Therefore, the present paper generalizes a study presented in Sec. Ill of HI, which 
is related to binary systems. It is worth stressing that the approach considered here 
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has been performed in a completely analytical way. The validity of this approach is 
supported by finding that the power lost in gravitational waves by the system agrees with 
the well-known result obtained by Peters and Mathews Pi using a different approach. 
In the next sections, the case of a binary system composed of two point particles of 
different masses moving in circular orbits in a Minkowski background is studied. In 
Sec. II, we present a brief review of the formalism, the notation used and the Einstein 
field equations in the characteristic formulation on null cones oriented to the future. In 
Sec. Ill the specific solutions for the field equations are constructed by expanding the 
metric variables in terms of the s Zi m spin-weighted spherical harmonics. A discussion 
on the boundary conditions satisfied by the coefficients in these expansions for arbitrary 
values of the l and m modes is provided. In Sec. IV the gravitational radiation emitted 
by the binary system is computed using the Bondi News function. In particular, the 
contribution to the gravitational radiation by the (7) odd modes due to the asymmetry 
(different masses) of the system is considered, and it is shown how these terms vanish 
if the system is symmetric (equal masses). Finally, in Sec. V other issues for further 
study are presented. 

2. Formalism 

The characteristic formulation based on future-oriented null cones, when the eth 
formalism is employed, is well known and has been used to explore a wide variety 
of problems [9], [13l 1~1. 12], [16) 20]. However, a brief review of some of its most important 
aspects is necessary to present the notation and conventions used here. We use the 
geometrised unit system, i.e. G = c = 1. The Greek indices run from 1 to 4, whereas 
the capital latin letters run from 3 to 4, labelling the angular indices. The coordinates 
used are = (u,r,x A ), where the retarded time u labels each null cone upon which 
the spacetime is foliated, r is the luminosity distance and x A represents the angular 
coordinates. Here x A can denote (#,0) or ( q,p ) for the usual spherical coordinates or 
for the stereographic two patch representation, respectively. In these coordinates the 
Bondi-Sachs metric takes the form 

ds 2 = _ ^2/3 (i + V) - r 2 h AB U A U B ) du 2 - 2 e^dudr 

— 2 r 2 h AB U B dx A du + r 2 h AB dx A dx B , (1) 

where w and the redshift /3 are related to the square of the ADM lapse function jlj. The 
shift vector between two successive null cones is labelled by U A . The metric associated 
with the angular manifold (unit sphere) is represented by h AB ( q A B )• In addition, the 
Bondi gauge is imposed, i.e. det(h AB ) = det(q AB ) = 1. Further it is required that 
h AB h BC = q AB q BC = 5ff. 

The metric of the unit sphere is decomposed as q AB = q( A q B ), where the overlinc 
indicates complex conjugation, the parenthesis symmetrization and the dyad q A are 
the null tangent vectors along the coordinate lines associated with the angular charts 
used to make the finite covering of the sphere. These vectors satisfy the relations: 
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qAQ A = 0 and qAq A = 2. In stereographic coordinates the q A vectors take the form 

\ / Q' 


/=(l + < 0(^ + ^ 4 >, £ — tan 


( = q + ip 


(see e.g. 12 El ED ED whereas in spherical coordinates (see e.g. panada]), 

q A = 5 a 3 + i esc 95 a 4 . 


All tensor quantities, such as the metric or the Ricci tensors, are contracting with the 
dyad producing spin-weighted scalars. The metric hAB for the angular manifold can be 
represented by J = q A q B hAB/2 and K = q A q B hAB/2 , and U A by U = q A UA■ Also, the 
Bondi gauge constrains K to satisfy K 2 = 1 + ././. Thus hAB depends only on J. It is 
worth stressing that the case J = 0 returns the symmetrical unit sphere metric. 
Associated with the projections of the covariant differentiation referred to qAB onto 
the q A vectors, the differential operators, labelled 5 and 9, are constructed. These 
operators can raise or lower the spin-weight of any spin-weighted scalar function. The 
spin-weighted functions are given by 

m n 

i=l j =1 

whereas 9 is defined 

m n 

s =? d ii a ^ n ^ A ^ Bi '" Bm M- An \ D > 

i= 1 j =1 

and 9 is defined 

m n 

s .*=. 

i=1 3= 1 

where the vertical line in the indices indicates covariant differentiation associated with 
qAB, and the symbols A s . and \ Ai can take the values q Bi or q B . and q Ai or q A \ 
respectively. Specifically, these operators take the following form 

9 a V = q D s^,d + sQ , 9 S T = q D S T D - sU , 

where = —q A q B q B \A/2- For a complete review of these operators and their properties, 


The Einstein equations can be written as 

Eyu, = Rftu - 8?r (T„ u - g^T/2) = 0 , 
and in this formalism, they are decomposed as 


E 22 — 0 , 

E 2A q A = 0, E AB h AB = 0, 

(2a) 

E A Bq A q B 

= 0 , 

(2b) 

Eu — 0 , 

O 

0 

CN 

&q 

(2c) 


which corresponds to hypersurface, evolution and constraint equations, respectively 

(HUEUmUTUIH]. 
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In the linear regime, the following set of equations corresponding to a perturbation in 
the Minkowski background, as originally obtained by Winicour [23] and later by Bishop 
na, is given by 


4/3, 


Sn Tq 2 = 


8vr T 2A q A = 


9 7, 


2 9/3 (r 4 U, r ) 

+ —- + 


8vr ( h AB T AB - r 2 T ) = -299/3 + 


r 2 r 2 ’ 

2 7+3 2 j (ryse + glT)) 


2 2r 2 

+ 4/3 — 2w >r , 

87rT AB q A q B = -29 2 ^ + (r 2 91/) r - (r 2 7, r ) + 2r (r J) ur , 


87T 

87T 


T \ 99 w 99/3 (3tf + 9*7) ™ ^ 2/3,, 

+ 7 11 I — „ o H-n---1-— + 


2 ) 2 r 3 r 2 2 r 2 2r r 


T 


+ r 12 = 


99,9 (r 2 (917 + 917)) r w 


4^2 


+ 


2r ’ 


8nT lA q = 


A 9 J u 9 2 f/ ,9917,1 ^9w\ „„ , {r A U) 


2 4 

r 2 U , 


+ 


4 + 2(t) “ 8 + + 2r 2 


' ,w 


+ u 


(3«) 

(36) 

(3c) 

(3d) 

(3c) 

(3/) 

(3s) 


3. Binary system 


It is worth stressing that one of our aims is to study the well-known problem of a circular 
binary system of point particles with different masses. We show that the Peters and 
Mathews result for the power radiated in gravitational waves (see Appendix A) can 


be obtained by using the characteristic formulation and the News function. Also, the 
present paper generalizes previous results d applying to particles with equal masses. 
The particles held together by their mutual gravitational interaction are far enough 
from each other such that to first order, the interaction between them can be considered 
essentially Newtonian. This assumption is valid if one considers the weak held 
approximation, in which the Bondi-Sachs metric in stereographic null coordinates 
reduces to, 


w 


ds 2 = — 1 —--2/3) du — 2(1 + 2/3)dudr — 2r 


2(U + U)^' n »i(U 


U) 


! + | C |2 d( l du ~ 2r "i + \ c ]2 d P du 


2 r 


(2 + 7 + 7) 

(1 + ICI 2 ) 2 


dq 2 — 4ir 2 


M 

(1 + ICI 2 ) 2 


dqdp — 2 r 


(—2 + 7 + 7) 

(1 + ICI 2 ) 2 


dp 2 


(4) 


Note that gu ~ —1 + 2$, where $ = /3 + w/(2r) represents the Newtonian potential, 
as usual in this kind of approximation. 

We take these two particles in a Minkowski background, analogously to Peters and 
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Mathews [12] and Bishop et. al. in id- Such a system allows one to explore in full 
detail the boundary conditions across the hypersurfaces generated by their orbits shown 
in figure HJ 



Figure 1. Binary system with the world tubes of each orbit extended along the 
direction of the retarded time, separating the spacetime into three regions. 


The density that describes the binary system is given by 


P = 


5(9 - n/2) 


Mi5(r — ri)5(4> — vu) + M 2 5(r — r 2 )5(f — vu — it) ) , (5) 


where, r* (Mi) are the orbital radius (mass) of each particle and r\ < r 2 . 

The orbits of the masses generate world tubes, which are extended along the retarded 
time coordinate, allowing the separation of the spacetime into three empty and disjoint 
regions: inside, between and outside the matter distribution. 

In order to solve equations fl3ajl -13o) for the vacuum, the following expansion for the 
metric variables is used 


s / = 53»(/ {m e i H w )g i 0 Z {m , (6) 

l,m 

where m means Y^ m =-v s ^ lm are ^ ie spin-weighted spherical harmonics, which 
are defined in [25] as s Rim and are eigenfunctions of the [5, 9] operator and s f represents 
the functions /3, W, U, U, J, J. The substitution of equation ([6]) into (]2]) provides a system 
of ordinary differential equations for the coefficients in the above expansions. The 
families of solutions, for l = 2, satisfying this system of equations for the vacuum read 


fomir) 

J2m(r) 


U 2m (r) 


Dlf)2m ; 

2iDil32m _ Dij 2 m(vr \m\ - 1 )(ur \m\ + 1) 
vr \m\ 6r 3 

_ iD 2 j 2 m e 2ivr \ m \vr \m\ + if T> 3J2m (z/r |m| - 3i) 
8 u 5 r 3 \m\ 5 vr \m\ 


2D lf s 2m (ur 

\m 

| + 2i) D 1J2m (2 v 2 r 2 

m 

2 + 4 ivr 

\m\ 

1+3) 

vr 2 

m\ 

6 r 4 


(7a) 


(7b) 
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U >2 m(r) 


D2j2me 2iurlm K2ur 

\m\ 

+ 3i) iD 3J 2 m ( u 2 r 

2 ' 

m\ 

| 2 + 6 ) 

8 u 5 r 4 

\m\ 

5 

ur 2 

\m 



lQrD ip2m + 6rD 3J2 m(2 + iur \m\) 


^D2J2me 2iUr ^ 
4 u 5 r 2 Iml 5 


(7c) 


_ iDu2m{(7 + i)vr \m\ - i)(l + (1 + fur \m\) ^ 

where the constants of integration are represented by D n Fim\ here n indexes the constant 
and F corresponds to the Bondi metric function whose integration generated it. This set 
of families of solutions depends only on four constants, namely, Dip 2m , D 3J 2 m , DiJ 2 m and 
D 2 .j 2 m- This is so because the families of solutions for the coefficients /32m, U 2m and 
w 2 m , resulting from ffTal) - ff3dl) . are constrained by (IHejh (pT/|) . This fact is independent of l, 
thus the set of families of solutions for any l will have four degrees of freedom. A unique 
solution for the whole spacetime cannot be determined by only imposing regularity 
of the metric variables at the null cone vertices and at the null infinity. Therefore, 
additional boundary conditions must be imposed. In particular, this can be done by 
imposing boundary conditions on other hypersurfaces, such as in the case of the thin 
shells studied by Bishop [15] , in which the additional conditions are imposed across the 
world tubes generate by the shell itself. Once the above constants are determined, one 
readily obtains the metric functions /?, J, U, and w for the whole spacetime. 

As divergent solutions are not expected at the vertices of the null cones, regularity 
conditions at these points must be imposed for the metric. In order to do so, an 
expansion of the metric variables around r = 0 in power series of r is made and the 
divergent terms are disregarded. This procedure establishes relationships between the 
coefficients, leading to a family of solutions for the interior that depends only on one 
undetermined parameter, where in particular /3; m _(r) = 0. One obtains, for example, 
for / = 2, 


P2m-(r) 
Jim— (j ") 


U 2m -{r) 


U>2 m-(r) 


0 , (8a) 

— ^ 2J2m ~ (2 u 3 r 3 \m\ 3 — 3 iu 2 r 2 \m\ 2 e 2lur \ m \ — 3iv 2 r 2 \m\ 2 
24 u 5 r 3 \m\ 

+6 ur \m\ e 2iurlml + 3ie 2ji/r|m| - 3 i) , (86) 

iD 2 J 2 m- (r, 4 4 | |4 . c 2 2 I |2 i i 2 ivr\m\ 

- (2u r \m\ + bu r \m\ — biur \m e 1 1 

24^ 5 r 4 \m\ 

—12iur \m\ + 9e 2ivr|m| — 9) , (8c) 

- fffSL — (2iz/ 4 r 4 |m| 4 + 4 u 3 r 3 \m\ 3 — 6 iu 2 r 2 |m| 2 — bur \m\ 

4 u 5 r 2 \m\ 

—3ie 2ivr \ m \ + 3i) . (8 d) 


For the intermediate region, the same structure of the general solutions is maintained, 
for the case l = 2 given by fl 7 o.p - (1 7 cifl . That is so because there is no reason to discard 
any particular term, or to establish any relationship between the constants as occurs 
for the interior region. For the exterior region, since that is required regularity at the 
null infinity, the coefficient of the exponential factor (exp(2iur\m\)) must be null. This 
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means that all constants of the form D 2 ji m +, with / = 2, 3, - - must be zero. Therefore, 
the number of degrees of freedom for the exterior family of solutions is reduced in one 
parameter. Thus, a family of solutions to describe the whole spacetime for the field 
equations (J3aj)-(J3^), with eight parameters to be determined is obtained. Now, in order 
to fix these eight constants, it is necessary to impose additional boundary conditions, 
in particular across the time-like world tubes generated by the orbits of the particles. 
These boundary conditions across the world tubes, i.e. when r — i — 1, 2 come from 
imposing discontinuities on the metric coefficients, 

l9u] ri = 0 > [ 912 ] n = Ag 12 \ r . , \giA] ri = 0 , [ 022 ] rj = 0 , 

[ 92A] ri = 0 , [g 3 „] r . = 0 , [g^] r . = 0 , (9) 

and on their first derivatives, 

= A 9^ , !*,” = i,"’ 4 , (10) 

where the brackets mean [f{r)\ r . = f(r)\ r . +e — f(r) | e . 

From the linearised Bondi-Sachs metric (ED) . and from the two sets of jump conditions 
(JH1) and PD , the coefficients /3; m , Ji mi Ui m and wi m are restricted to satisfy 

[^im(Fj')] AWji m ) \filmij'j)\ A/3ji m i 

[Jlm(rj)\ = 0 , [UlmiTj)] = 0 , (H) 

and for their first derivatives 


K m( r j)} = , [Plm( r j)} = A Pjlm , 

[J'im( r j)} = AJ' Im , [ULirs)] = A U’ lm , (12) 

where j = 1, 2, and A w j[mi A/3 jZm , A w' jlm , AA J' jlm and AU' jlm are functions to be 
determined. 

Solving equations (1111) and 03, simultaneously for both world tubes, the boundary 
conditions are explicitly obtained. We find that 


A(3jim bji m , Awji m 2 rjbji m , 


(13a) 


where bji m are constants. Note that this last fact implies that A f3'j lm = 0. We determine 
that the jumps for the first derivative of the J lm and Ui m functions are given by 


AT/ _ 8 v 2 rjbj lm H_ 

jlm (l- 1)1(1 + !)(/ + 2) ’ 


A Ujim — 2bn m 


4iv\m\ 


(13 b) 

1(1 + l)r ( ,' ' (13£,) 

The boundary conditions (11361) and (113 d) fix all parameters of the families of solutions, 
providing the specific solutions for the coefficients j3i m , Ji m , Ui m and wi m . These 
coefficients can be written as 


fim(r ) = fum(r) (1 - 0(r - r x )) + f 2lm (r) (0(r - n) - 0(r - r 2 )) 
+ f 3 im(r)<~>(r - r 2 ) , 


( 14 ) 
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where fi m represents /3/ m , J; m , Ui m and Wi m , with the first subscript on the right hand 
side terms indicating the interior (1), the middle (2) and the exterior (3) solutions; and 
@(r) is the Heaviside function, namely 


0(r) 


0 r < 0 

1 r > 0 


(15) 


These solutions depend explicitly on two specific parameters, namely bji m , with j = 1,2, 
which are related to the density of matter. The specific form of these relationships is 
obtained by just integrating the first held equation (l3al) across each world tube. As a 
result one obtains 


hjlm 2 TtrjPjim (l T te) , (16) 

where pji m are given by 

Pjim = - [ d(uu) [ dZl [ dr o Z^e-^p , (17) 

^ Js Jn Jij 

in which S = [0, 27t), v } is the physical velocity of the particle j in the space, and Ij is 
an interval e around r : that is given by Ij = (rj — e/2, r 3 + e/2), with e > 0. 

Before proceeding, it is worth noticing that the above procedure is a generalization of 
Section 3 of the paper by Bishop et. al. D3, in which the binary components have 
equal masses. In particular, the boundary conditions are also generalized since in the 
present case there exist two independent world tubes. Interestingly our solution is fully 
analytical. 

In order to include the null infinity, which is reached when r tends to infinity, a radial 
compactihed coordinate s is defined as follows 
r 

r + R 0 

where Ro is a compactihcation parameter. Thus, 0 < s < 1, where s = 0 and s = 1 
corresponds to the null cone vertices and the null infinity, respectively. 

Figure [2] shows some of the coefficients of the expansion of the metric variables in terms 
of the compactihed coordinate s for l = m = 2. 


M,= 1 A4=£ r 1 = ! r 2 = | fib =2 v= df M, = 1 W = | n = | r 2 = f flb =2 v= £- 
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M|=1 n = ! /2=| fio = 2 „ = a£ 


M=1 A4=4 /-I = | f2 = | Ro = 2 u = i/l 


•sr 3 


-i 


4 


5 


0 



0.0 0.2 0.4 0.6 0.8 1.0 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


S 


S 


(d) 


Figure 2. Real part of some components of the metric functions (l = m = 2) versus 
the compactified coordinate s (see the text) for a binary system with Mi = 1/2, 
M 2 = 1. The frequency of rotation is computed by means of Kepler’s third law. Here 
r\ and T 2 are referred to the center of mass of the system. 


Here Mi = 1/2, M 2 — 1, Ro — 2, and the radius of each orbit is referred to the center 
of mass of the system, namely 



where /1 is the reduced mass of the system and d 0 is the distance between the masses. 
The frequency of rotation v is computed by means of Kepler’s third law, i.e., 



(19) 


It is worth noting that the jumps in and Wi m functions are present at exactly rq 
and r 2 , whereas for Ji m and Ui m only their first derivatives present discontinuities, in 
agreement with the boundary conditions (1136|) and (113 eft . 

To illustrate the behaviour of (3, J, U and w we present them in figure [3] as a. function of 
s and 0 for a particular value of the retarded time u. These functions are constructed 
by using equation (J6]), and the solutions for the coefficients for each l and m. In this 
case, we use l < 8. As expected, the metric functions 0 and w and the first derivatives 
of J and U show jumps at (r, 6, 0) = (rq,7r/2, uu) and (r, 0, 0) = (r 2 , 7 t/2, vu — ti), which 
are just the positions of the masses, in agreement with the boundary conditions initially 
imposed. 

Note that since the first field equation for the vacuum 0 >r = 0 implies that 0/ m are 
constants along r, as sketched in figure [21 and that 0 is a gauge term for the gravitational 
potential, then, $ can be redefined as $ = w/(2r). These facts make the choice of the 
frequency v as obeying Kepler’s third law completely consistent and natural. 
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(c) 


(d) 


Figure 3. Snapshots of the metric variables as seen from the equatorial plane 
(9 = 7t/ 2), as a function of s and <p for u = 7r/2. Here Mi = 16, M 2 = 4/3, r\ = 1/13, 
r 2 = 12/13, Rq = 1/2 and v = 2-^/l3/3. (a) {3{s,<j>), (b) J(s,<j)), (c) U(s,<j)) and (d) 
W{s,4>) = w{s,4>){ 1 - s 2 )/(s 2 Rl). 


4. Gravitational Radiation from the source 


Now, we proceed with the calculation of the power emitted by the binary system via 
gravitational wave emission. We show that the approach presented here is consistent 
with the well-known result obtained by Peters and Mathews [19] . 

Following Bishop [15] , the Bondi News function in the weak field approximation is given 

by, 


u 


lim 

r —>00 


r 2 J, 


+ 


d 2 o; 



2 


2 
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Substituting the metric variables given in (15|) . one obtains the News function for l > 2 
and —l<m<l, namely 

U = lim V s [ + ‘1 + g 2 oZlm . ( 20 ) 


l,m 


Now, substituting the coefficients of the metric variables for the exterior region, one 
obtains 


a r= - 


in 3 2 * 5*21 .. /2 3 c in 4 2 S 31 

a /6 V 3 a/30(z/ — 3i) 


8 i 


2 ,.4 C 
— Z/ 2 *> 


2*^32 


9* 


AZ2 4 2^33 


iv 5 9 S', 


2041 


2z/ — 3i 
%i\[\v h 2*5*42 


^ - i 3 a/T0 (z/ 2 - 7iv - 14) 3 (2z / 2 - 7iv - 7) 

81 iv 5 26*43 


2567/t^ 5 2 6 , 


2*^44 


y/io ( 9 z 2 2 - 21 iz 2 - 14) 3 ( 8 z / 2 - 14iz/ - 7) 

2 S'; 


2*J51 


+ 


16 \/ yj ^ 6 2*S*52 


a/2!o (iv 3 + 12 z 2 2 - 54 iv - 90) 4 iv 3 + 24z/ 2 - 54iz/ - 45 


27a / 2*S*53 


1024,/ Z 2 6 2 S'; 


2*^54 


3 iv 3 + 12z / 2 - 18iz2 - 10 32iv 3 + 96z/ 2 - 108 iv - 45 


625 


5 . .6 c 
— V 2 *> 


2*^55 


+ • • • , 


( 21 ) 


( 22 ) 


25 iv 3 + 60 z / 2 — 54 iz 2 — 18 
where we define the spin 2 quantity 2 S) m as 

_ (K(ZW+e iH ™) S 2 o7z m + ^(Duz-m+ZH™) 6 2 0 7z _ m ) 

2 y/(l-m + l)(l + 2) 

Since the binary system is confined to a plane, then a natural choice to simplify the 
problem of expressing the News function, without loss of generality, is to put the masses 
motion in the equatorial plane 6 = 7 t/ 2 . This means symmetry of reflection for the 
density of matter and, consequently, for the spacetime. Thus, this choice restricts the 
components of the density, obtained from (fTTlh to have the following form 

M^rfS (r — r 2 ) + Mir\5 (r — ri) 


Plm 


Plr. 


Plr 


2 ir ,2 


rrr. 


1' 2 


Mir Id (r — ?•() — M 2 r 2 <5 (r — r 2 ) 


2^2 


rrr, 


l : m even 


l,m odd 


(23) 


1 ' 2 


where pi m are numerical constants. Therefore, for binaries of different masses, the News 
function fOTj) is simplified to 


M — — / ‘fv 3 2 * 5*22 - 


iv 4 2*S*31 

a/30(z2 - 3i) 


9 i 


3 , .4 c 
V 2*J; 


2*J33 


8A/|z2 5 2*51 


2*->42 


V — l 


3 (2v 2 - 7 iv - 7) 


2567/V 2 S, 


2*^44 


+ 


Z2 6 2 S'; 


2*- , 51 


3 ( 8 z / 2 - 14w - 7) ^210 (in 3 + 12z / 2 - 54 iv - 90) 
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27A / oS, 


2^53 


625\/J^z/ 6 2^55 


3 iu 3 + 12z/2 - 18 iu - 10 25 iu 3 + 60z/ 2 - 54iv - 18 


+ 


• (24) 


When the explicit solutions are used, the News functions for the binary system take the 
form 

' 1 


[2n 

M = 8 y yr - 2-^22 (M2] + A 4 


22) v" + —iJ — 2-^31 (AI31 — M32) V 4 


I3tt 


9 i\ — 2L33 {-M -31 — -M32) Z2 4 + — V / 27T 2-^42 (/W41 + /W42) V° 
7 63 


~TT \ ~TT 2-^44 (/W 41 + .M 42 ) Z' 5 — —i\! -yrr 2-^51 (2^51 — 2 W 52 ) Z2 6 

9 V 7 180 V 164 

27. [3tt . 6 625. [bi r . . 6 

-— 2-^53 (A 4 5 i — M.52) V + — 7 T l \ Xy 2-^55 (.M51 — .M52) v 


40 V 11 


24 V 33 


where 


+ 1 ) , 


(25) 

(26) 
(27) 

Note that, as consequence of (1231) . for Mi = M 2 = M 0 the terms with l odd disappear 
from the News function (1251) . Thus, as expected, one obtains immediately 


and 2 Lim are defined as 

2 L lm = ( 2 Z l _JR(e‘H-) - \™) 2 Z lm ) 


/V — 16 \j —z/ 3 Morg (Vq + l) 2-^22 + ~V / 27rz y5 114orQ (l/ 0 2 + l) 2-^42 


~\j~~^ uJM o r o (K) 2 + l) 2-^44 + X485 u7 M°r° ^ + 1 ) 2 7/62 


8192 nr 


495 


y^M„r$(A? + l) 2 L« + lfl V /^^M 0 r»(V? + l) 2 L« 

(28) 


where Vo is the physical velocity of the identical masses, which is obviously tangent to 
the circular orbit. 

The energy lost by the system dE/du is related to the News function, via 


if = L f dn Mjr , 

du 4vr J n 


(29) 


which results for Mi f M 2 in 

= ^6 (M21 + _A/f22)^ + ^ (A 4 3 i - M 32 ) 2 

57376 io/,2 a2 4010276 12 2 

+ 1969 " + + 755925"" (Msl ~ M * ] + "'' (30) 

Notice that the hrst term on the right side of the above equation is identical to the 
power lost obtained by Peters and Mathews [T9] for circular orbits, whereas the other 
terms stand for the octupole, hexadecapole, etc. 
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5. Summary and Conclusions 

In this study we generalize a previous work by Bishop et. al. mi, concerning equal 
mass binary systems in a Minkowski background. Here we consider the case in which 
the components of the binary systems have different masses, although still in circular 
orbits. 

We show that, instead of two regions, as in the case of binaries with equal components, 
the spacetime needs now to be separated in three regions: interior, between and outside 
the two world tubes. As a result, the match conditions satisfied by the coefficients in the 
spin-weighted spherical harmonics expansion, across the two hypersurfaces generated by 
the circular orbits of these two (different) masses were generalized. Also, the procedure 
developed here allows one to perform calculations for arbitrary values of the l and m 


modes. 


It is worth stressing that one of the interesting aspects of the present study has do with 
the development of a procedure that can be applied in problems in which multiple layers 
are present. 

We also calculate the energy lost by the emission of gravitational waves by means of the 
Bondi News function. Again, we do that for arbitrary multipoles, in other words, for 
different values of the l and m modes. The interesting point here is that for different 
masses the emission of gravitational radiation occurs for all values (multipoles) of l > 2; 
for the particular case of binary systems with equal components, the multipole terms 
for odd values of l vanish. 

Last, but not least, the present study can be used as a method to extract gravitational 
radiation of binary systems during their periodic phase, when one considers that the 
system can be surrounded by a boundary, T, very far from the system, where a linear 
approximation can be applied. 
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Appendix A. 

It is worth recalling the main steps given by Peters and Mathews (19j to obtain the 
well-known and widely used equation for the power radiated by two point masses in a 
Keplerian orbit. 

In the weak field limit of the Einstein field equations, i.e., when the metric can be 
written as a perturbation h^ of the Minkowski metric rj^ u , namely 



(A.l) 
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the power emitted by any discrete mass distribution in the limit of low velocities, as 
shown in [19], is given by 

P = 5 (pijQij — 3 ’ (A-2) 

where the dots indicate derivative with respect to the retarded time u and, 

Qij ^ ' TflaXaiXaj i (A.3) 

a 

where a labels each particle of the system and x a i is the projection of the position vector 
of each mass along the x and y axes. Particularly, for a point particle binary system of 
different masses in circular orbits, when the Lorentz factor is considered as 7 = 1, one 
can write 

x a i = r a cos {vu - 7 r< 5 a2 ) , x a2 = r a sin {vu - n 5 a2 ) , a = 1 , 2 , (A.4) 

as shown in figure IA1I 



Figure Al. Binary system as viewed from the top. The coordinates x a i of the 
particles are indicated, as well as the angle vu with respect to the x axis. 


Qij 


Here, r a is given by f[T 8 ]h v by (fT9]h and the Kronecker delta distinguishes each particle. 
Then, the components of reads 

/adl cos 2 {vu) fidl sin(Vw) cos (vu) 
j-jjd/Q sin(z/n) cos (vu) /idg sin 2 (vu) 

thus, 

- f Av^jidf sin( 2 z/u) —4z/ 3 /idQ cos( 2 z/u) 
lj \—Av z ndfcos{2vu) —4j/ 3 /idQ sin( 2 z/u) 

Finally, substituting the above equation in (IA.2jh one obtains 


(A.5) 


(A. 6 ) 


„ 32 2 6 4 32mi 2 m 2 2 (mi + m 2 ) 

P = T '“ /<i » = - 54 - 

where Kepler’s third law is used in the last equality. 


(A.7) 
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